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Abstract 

We consider some lattices and look at discrete Laplacians on these 
lattices. In particular, we look at solutions of the equation 

A(1)0 = A(2)Z, 

where A(l) and A (2) denote two such Laplacians on the same lattice. 
We show that, in one dimension, when A(i) i = 1, 2 denote 

A(l)^ = ^(i + 1) - - 1) - 2^(z) 

and 

A(2)Z ^ Z{i + 2) - Z{i - 2) - 2Z{i) 
this equation has a simple solution 

(j>{{) ^ Z{i + 1) + Z{i-l) + 2Z{i). 

We show that in two dimensions, when the system is considered on a 
hexagonal (honeycomb) lattice, we have a similar relation. This is also 
true in three dimensions when we have a very special lattice (tetra- 
hedral with points inside). We also briefly discuss how this relation 
generalises when we consider other lattices. 



* email: w.j.zakrzewski@durham.ac.uk 
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1 Introduction 



Recently, a Prohlich Hamitonian was studied on two-dimensional, discrete, 
quadratic [1, 2, 3], and hexagonal lattices [4]. The resultant equations were 
rather complicated to solve but when one restricted oneself to looking at sta- 
tionary fields, the equations simplified somewhat and, in the case of hexago- 
nal lattice, one could decouple them and reduce them to the localised discrete 
nonlinear Schrodinger equation. 

The reason for this decoupling lies in the observation that, on a hexagonal 
lattice, the equation 

A{l)(f) = XA{2)Z (1) 

has a simple localised solution for Z as function of (p. Here A(z). i = 1,2 
refer to the two simplest discrete Laplacians on a hexagonal lattice - i.e. 
those involving 4 and 7 lattice points (more details will be given later). 

A similar situation holds in one dimension and, in fact, was used by 
Davydov [5] in his observation that the interaction of an appropriate 
Schrodinger field, such as a field describing amide I- vibration in biopoly- 
mers, with the distortions of the underlying lattice, results in the creation of 
a localised state which has, since then, been refered to as Davydov's soliton. 

In this paper, we look in detail at the eqution (1) in the case of various 
lattices. We find that a regular lattice in one dimension, a hexagonal lattice 
in two dimensions and a tetrahedal lattice (with points inside) in three 
dimensions are somewhat special as it only in their cases the system possesses 
simple localised solutions. 

We discuss the modifications that are required to have localised solutions 
for other (more complicated) lattices. 

2 One dimension 

Consider a regular lattice as shown in Fig. 1. 



i-l I i+1 

Figure 1: One dimensional lattice. 
For simplicity let us take the lattice spacing to be given by a = 1. Then 
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define A(l)P(i) as 

A(l)P(i) = P{i + 1) + P{i - 1) - 2P(i). (2) 

and 

A(2)P(z) = P(z + 2) + P{i - 2) - 2P(z) (3) 

i.e. the same expression as above but with the displacement by two lattice 
units. 

Then consider the equation: 

Aimi) = A(2)Z(i), (4) 
which we want to solve for ^(i) in terms of Z{k). 

2.1 Continuum limit 

Note that in the continuum limit our equation (4) reduces to 

— — = 4 

which clearly has a solution 

(t){x) = AZ{x) 

plus, of course a linear and a constant piece. 

2.2 Discrete case 

In this case it is easy to see that a solution is given by 

(t>{i) = Z{i + 1) + Z{i - 1) + 2Z{i). (7) 

Of course, to this we can also add Ai + 6, which are the lattice equivalents 
of the linear and constant pieces of the continuum case. 

Note that our solution (7) has the correct continuum limit. Note also 
that its structure is very similar to the structure of A0(z) - except that, this 
time, all terms carry the + sign. 



3 



(5) 
(6) 





















(i.j+1) 








(i-l, 


) 


(ij) 












(iJ-1) 































Figure 2: Two-dimensional square lattice. 

3 Two dimensional Cases 
3.1 Squeire lattice 

When we have a square lattice (see fig. 2) we can consider various Laplacians. 
One of them, which we shall call A(l), involves the nearest lattice points 
and is defined as 

A(l)0(z,j) = cl){i + lj) + 4>{i-lj) + cl)(ij + l) + cj)(ij-l)-mi,j). (8) 

The second 'obvious' Laplacian can be defined as the one above with shifts 
by two lattice points ie 

A(2)</>(i,i) = (^{i + 2,j) + ,pii-2,j)+^{i,j+2) + <Pii,j-2)-mi,j), (9) 

or we could also 'involve' the corners (ie terms corresponding to i=b 1, j ± 1). 
Thus we could define A(2')^(i,j) as 

A(20</>(i, j) = 0(i + 2, j) + ^{i - 2,j) + cp{i,j + 2) + <p{i,j - 2) 

+a (0(z + (l>{i + 1, j - 1) + 4>{i - 1, j + 1) (10) 

+,^(i-l,i-l))-4(l + a)<^(z,j), 

for any value of a. 
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The work mentioned in [1, 2, 3] required a solution of 

A{l)P{ij) = fiA{2)Z{ij), (11) 

where A(2) is given by the expression as in (13). 

Unfortunately, this equation has no simple local solution. The problems 
is with the 'corners'; the second Laplacian should involve expressions at 
i±l,j + ±l. 

Thus it is possible to solve 

A(l)P(ij) = tiA{2')Z{ij) (12) 
for an appropriate choice of a. To see this, in analogy with (7) take 

P{ij) = X{Z{i-l,j) + Zii + l,j) + Z{i,j-l) + Z{i,j + l) + AZ{i,j)). (13) 
Then 

A{l)P{i,j) = X{Z{i + l,j) + Zii-2,j) + Z{i,j-2) + Z{i,j + 2) (14) 

+2Z{i + 1, j + 1) + 2Z{i - l,i + 1, +2Z{i + 1, j - 1) + 2Z{i -1,3- 1)) , 
which is clearly a solution when a = 2 in (10). 

3.2 hexagonal (honeycomb) lattice 

Next consider a hexagonal lattice as shown in Fig.3. 




Figure 3: Hexagonal (honeycomb) lattice. 
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For A(l) take the 4-point Laplacian ( involving and 3 'nearest' points); 
ie take 

A(1)P(0) = P(a) + P{b) + P(c) - 3P(0) (15) 

For A (2) take the Laplacian involving the 'next to the nearest' points - 
ie. the 7-point Laplacian defined by 

A(2)Z(0) = Z{1) + Z{2) + Z{3) + Z(4) + Z(5) + Z{6) - 6Z(0). (16) 

Then consider 

A(1)P(0) = /xA(2)Z(0). (17) 

In fact, this is the equation which arose in the study reported in [4]. There 
it was shown that this equation has a simple local solution and that solution 
is given by 

P(0) = ^{Z{a) + Z{b) + Z{c) + 3Z(0)). (18) 

Thus, in this sense, the hexagonal lattice resembles the most, the one 
dimensional case; the double lattice spacing in the one diemensional case is 
replaced by the double the angle between the directions on the lattice. 

3.3 trianguleir lattice 

Next we consider a triangular lattice as shown in Fig. 4. 

This time the obvious A(l) involves taking taking the 7-point Laplacian 

A(1)P(0) = P(l) + P(2) + P(3) + P(4) + P(5) + P(6) - 6P(0). (19) 

For the 'second' Laplacian we have more choice. We can take 

A(2)Z(0) = Z{a) + Z{h) + Z{c) + Z{d) + Z{e) + Z{f) - 6Z(0). (20) 

or we could add to it 

+ A{ Z{u) + Z{v) + Z{w) + Z{x) + Z(y) + Z{z) - 6Z(0)). (21) 

However, in general, it is difficult to solve the equation 

A(1)P(0) = A(2)Z(0) (22) 

Only when A = 2, i.e. for A(2') given by (24) and (25) with ^ = 2, we have 
succeeded in finding a local solution. 
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Figure 4: The triangular lattice 



This solution is given by 

P(0) = Z{1) + Z{2) + Z(3) + Z(4) + Z(5) + Z(6) + 4Z(0). (23) 
If guided by the experience gained from the previous cases we take 

P(0) = Z{1) + Z{2) + Z(3) + Z{A) + Z{b) + Z{<o) + 6Z(0). (24) 
we find that the Z field satisfies 
Z{a)+Z{b)+Z{c)+Z{d)+Z{e)+Z{f)+2 [Z{1) + Z{2) + Z(3) + Z(4) (25) 

+Z(5) + Z(6) + Z(u) + Z{v) + Z{w) + Z{x) + Z{y) + Z{z)] - 36Z(0) = 0. 

We see that in this case we have 19-point Laplacian - involving both the 
new and the original lattice points. 

One can consider other, more complicated, cases but it is clear that the 
hexagonal lattice is very special as only for it, as in the one-dimensional 
case for a regular lattice, we have a simple local solution of the problem 
(17) with Laplacians involving the 'nearest' and 'next-to-the-nearest' lattice 
points only. 
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4 Three Dimensions 



The discussion of the previous sections generahses very easily to three dimen- 
sions. First of all, for a regular cube-like lattice, for which in 2 dimensions 
we had problems with 'corner' terms, the situation is similar except that 
now we have more such problems. 
Thus if we take 

A{l)P{i,j,k) = P{i+l,j,k)+P{i-l,j,k)+P{i,j-l,k)+P{i,j+l,k) (26) 

+P{i,j, k-l) + P{i,j,k + l)-QP{i,j.k) 
and then consider 

P(i, J, k) = Z{i + l,j,k)+ Z{i -l,j,k) + Z{i,j -hk)+ Z{i,j + 1, k) (27) 

+Z{i,j, k-l) + Z{i, j, k + l) + 6Z{i, j.k) 

we find that Z{i,j, k) satisfies a rather complicated expression involving the 
fields at all points (i, j, k) which correspond to the shift of only one of these 
indices by ±2. There are 6 such terms and they all come with the coefHcient 
= 1. In addition we also have to add all terms which involve the shift of two 
of the indices {i,j,k) by ±1. There are 12 such terms and we have to add 
them all with the coefficient = 2. And then finally we subtract — 24Z(0, 0, 0). 
A little thought shows that this is the obvious generalisation of the square 
lattice which we discussed in the previous section. 

However, we have a regular lattice which is a 'natural' 3-dimcnsional 
generalisation of the hexagonal lattice in 2 dimensions. Such lattice can be 
constructed in the following way: 

Take a regular cube lattice and delete every other lattice point as in- 
dicated in Fig. 5. Then place a similar lattice at the points which are the 
centres of the cubes of the previous lattice. This way we have a body centred 
lattice with many points missing. 

Consider then the points of the original cube lattice as points of type A 
and those of the additional lattice as points of type B. Then the nearest 
neighbours of each point A are 4 lattice points of type B and vice-versa. 
The next to the nearest neigbours (ntn), however, come from the original 
lattice; ie for a point type A - are points of type A, and for the point B are 
points of type B. Each point A has 4 nearest neighbours and 12 ntn points. 
And, by construction, this is the case for both types of lattice points. 

Clearly we have two obvious Laplacians that we can consider - those 
involving the nearest neigbours and those involving only the ntns. 
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Figure 5: The three-dimensional lattice. Points A shown. In the centre of 
each other cube there is a point B. Points B, togetherm form another lattice 
- like the one shown but displaced by ^ lattice spacing in all directions. 

Let us consider, for simplicity, the point B located at (0,0,0) inside a 
cube of size 2a, with points, of type A located at (1,-1,1), (-1,-1, -1) and 
(1,1,-1) and (-1,1,1,) (we are using the convention of x - horizontally to the 
right, z vertically up and y away from the observer). Then we can define 

A(1)P(0, 0, 0) = P(l, -1, 1) + P(-l, -1, -1) + P(l, 1, -1) (28) 

+P(-1,1,1)-4P(0,0,0). 

In the continuum limit this expression, clearly, reduces to the Laplacian 
of P ie {dl + d^ + dl)P. 

For the second Laplacian we take ntn points. Notice that they correspond 
to points ±2 and such that one of the coordinates of x, y and 2 = while 
the others take values ±2. Thus we define 

A(2)P(0, 0, 0) = P(-2, -2, 0) + P(2, -2, 0) + P(-2, 2, 0) + P(2, 2, 0) 

+P(-2, 0, -2) + P(-2, 0, 2) + P(2, 0, -2) + P(2, 0, 2) (29) 
+P(-2, -2, 0) + P(-2, 2, 0) + P(2, -2, 0) + P(2, 2, 0). 
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Now we consider our equation: 

A(1)P = aA(2)Z, (30) 

where P and Z are taken at the same lattice point, say (0, 0, 0). 

Guided by our experience from the lower dimensional cases we take 

P(0, 0, 0) = -1, 1) + Zi-1, -1, -1) + Z(l, 1, -1) (31) 

+Z(-1,1,1) +4Z(0,0,0)). 

Then, it is a matter of simple algebra to check that P does indeed satisfy 
(30) when /? = a. A little thought shows that our lattice is really tetrahedral 
in nature (with points inside it), hence it is a natural generalisation of the 
hexagonal lattice in 2 dimensions. 

5 Conclusions 

In this short note we have discussed various lattices and Laplacians defined 
on these lattices. Our main interest was the relation between various Lapla- 
cians and the existence or not of fully localiseed solutions of the equation 

A(1)P(0) = A(2)Z(0) (32) 

involving these Laplacians. 

We have found that, in one dimension, this equation a localised solution 
when the lattice is regular and A(l) and A(2) involve the Laplacians con- 
structed with the 'nearest' and 'next-to-the-nearest' lattice points. Moreover 
the points in each Laplacian are equidistant from the central point. 

These conditions are also required in higher dimensional lattices where 
they are much more stringent. Thus in two dimensions they require that 
the lattice is hexagonal (honeycomb) and then A(l) and A(2) involve 1+3 
point and 1+6 point Laplacians. For other lattices the Laplacians mix points 
of different distance from the central point and we have not succeeded in 
finding a localised solution of (15). 

In three dimensions the relevant lattice is the 'body centred tetrahedral 
lattice'. In this case A(l) involves 1+4 lattice points and A(2) - 1+12 
points. 

Note that in 1 dimension - our Laplacians involve doubling of distance, 
in 2 dim - doubling the angle, and in 3 dim - doubling the spherical angles. 
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Noting the pattern of our results we expect that our results generalise to 
higher dimensions, where in D dimensions we expect A(l) to involve 1+{D+ 
1) points and A(2) 1 + D{D + 1) points. But lattices with such Laplacians 
do not appear to be physical so we have not attempted to construct them. 
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